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ABSTRACT 



The properties of an incompressible fluid in the mixing 
chamber of a constant area cylindrical ejector are analysed 
in a two-dimensional form. The flow field for velocity and 
temperature is assumed to be made up of regions in which 
uniform flow of the primary and secondary fluids exist, and 
which diminish with axial distance from the ejector entrance. 
The velocity and temperature distributions in the mixed 
region are assumed to have cosine shaped profiles. The com- 
patible solutions are given an axial distribution by assuming 
a parabolic spread rate for the secondary fluid jet boundary. 

Results are generated in the form of a non-dimensional 
velocity and pressure, a normalized temperature and a Momen- 
tum Factor, as functions of axial distance. Area ratios 
studied include 100:1, 9:1 and 2. 2£:l. Velocity ratios 
studied include 50:1, 10:1, 3sl and 1.5:1. The Fortran pro- 
grams employed to generate compatible solutions, titled 
EJECTMIX I and EJECTMIX II, aro included. 
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1. Introduction 

The majority of analyses of the physical properties of 
fluids in ejectors have used a one-dimensional approach in 
which uniform properties are assumed at the entrance to the 
mixing chamber which is of sufficient length for the pro- 
perties to again become uniform at the exit.-*- For example, 
consider a primary fluid flowing along the central axis of a 
cylindrical ejector and a secondary fluid flowing parallel to 
and completely surrounding it, as shown in Fig. 1. These 
fluids are separated by a circular boundary. At the entrance 
to the mixing chamber the boundary between the two fluids is 
removed and, at this initial cross-section, both the primary 
and the secondary fluids have uniform but different velo- 
cities. As the primary and secondary fluids flow axially 
through the mixing chamber, the fluids mix and the physical 
properties of each one of them are affected by the physical 

p 

properties of the adjacent fluid. 

It is the purpose of this investigation, by use of a 
high speed digital computer and Fortran Programming Language, 
to predict what the properties in the partially mixed, non- 

^•3 . Pai, Fluid Dynamics of Jets, D. Van Nostrand Go., 
Inc., New York, 19 5k 

^R. A. Smith, "Theory and Design of Simple Ejectors”, 
Some Aspects of Fluid Flow, Edward Arnold and Go., London, 
19E>1 
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;Fig,l Ejector Geometry and Initial Flow Parameters. 
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Fig. 2 Typical Velocity and Temperature Profiles 



uniform sections of the stream will be.^ Consideration of 
the fluid flow mechanisms, such as the shear stress within 
the fluids, leads to the concept of a profile shape for the 
properties at a given cross-section. A typical profile shape 
for velocity and temperature is shown in Pig. 2. 

The first problem considered is one in which the velo- 
cities of the primary and secondary fluids are different but 
their temperatures and densities are the same. The Fortran 
Program developed to solve this problem is called EJECTMIX I. 

The second problem, considered in the program EJECT- 
MIX II, studies the variations in the velocities and the 
temperatures of the two fluids. 



D. McCracken, A Guide to FORTRAN Programming, 
John Wiley and Sons Inc., New York, 1961 
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2. Ejector Analysis 

In order to analyse the behavior of the physical proper- 
ties in the mixing chamber of an ejector there are several 
possible analytic approaches. 

The first possible approach is a one-dimensional over- 
all analysis in which uniform properties are assumed at the 
entrance and uniform properties of the completely mixed 
fluids are assumed at the exit of the mixing chamber. With 
this analysis, the region in which mixing of the primary and 
secondary fluids takes place is not considered. The only re- 
gions where the physical properties are studied are at the 
entrance to the mixing chamber and after the point where uni- 
form flow exists. 

A second possible approach is that of an internal analy- 
sis which is based on the actual fluid mechanisms. This ap- 
proach requires the examination of an element of fluid and 
the forces which are acting upon it. The difficulty here 
lies in the inability to acquire a usable analytic expression 
for the shear stress acting upon the element of fluid, 

A third approach, the one used in this thesis, is the 
generation of compatible solutions for the physical proper- 
ties of the fluids at representative axial cross-sections by 
means of consideration of Conservation of Mass, Momentum, 
Bernoulli's Equation and the Energy Equation. In order to 
apply the above equations to the fluid in the mixing chamber, 
certain basic assumptions as to the nature of the flow in the 
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mixing chamber must be made. 

First of all, the flow in the mixing chamber is assumed 
to be divided into regions which are shown in Fig. 3a and 
Fig. 3c. These regions are defined by the central core 
boundary and the expanding jet boundary. There are similar 
velocity and temperature profiles which consist of uniform 
portions at the central axis and the periphery of the mixing 
chamber, which are connected by a cosine shaped profile. The 
assumed profile shapes are shown in Fig. 3b and Fig. 3<i» 
Asstiming profile shapes, negligible boundary effects, con- 
stant pressure across any given cross-section, and constant 
temperature until a mixing region is reached, the previously 
mentioned equations can then be used. 

The first equation considered is the Equation of Contin- 



uity 



m- f 3 p V dA - Constant 



For a typical profile shape, as shown in Fig. 3©» the 
Equation of Continuity becomes 



ra « 



>/ b f 

r a L 



Vprdr+/* V a ve+dV cos 



r-r. 



rb-r a 



-TT 



> S \ 

r dr+f Vs r dr] 



Integrating and evaluating the integral gives 

ra * T(p(vv r a 2 +Vave jr^-r^ -ij.dV +Vs |r g 2 -r b ‘ 

2 

Dividing through the entire equation by tr^Vgr^ , the 
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Typical Profile 
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in Region II 
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in Region III 



Pig. 3b Typical Profiles for Regions of Pig. 3a 
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Typical Profile 
in Region I 





in Region II * 



in Region Iil 



Fig«33 Typical Profiles for Regions of Pig. 3c 
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Temperature 
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Pig. 3© Typical Profiles with Equation Parameters 



equation takes on the non-dimensional form 
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T 



s r a 

o o 
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•p = -^ 

2 V 

S o 



r. V ave„ ^ r b« 2-r a^- 



a x , ave 
r 2 

a o 



a* 

r 2 
a 

© 



-ij-dVx 

V. 




The next consideration is that of Bernoulli’s Equation 
along a streamline in the uniform velocity portions of the 
profile. Bernoulli’s Equation states that 

2 

£ , V_ = Constant 
P + 2 



In non-dimensional form the equation becomes 

p , V 2 -Constant 
B f v s 0 v s 0 

The constant associated with Bernoulli’s Equation in the 
primary core is found by evaluating Bernoulli’s Equation in 
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the non-dimensional form along the axis at the entrance to 
the mixing chamber. Similarly, the constant associated with 
the secondary cone is determined at the entrance to the mix- 
ing chamber. These values will be used in conjunction with 
the Momentum Equation. 

Considering the same profile which was used in the 
Equation of Continuity and solving for the Momentum rate 



over the entire cross-section the Momentum Rate Equation be- 



Momentum rate 




comes 





Integrating and evaluating the integral 



2 2 

Momentum Rate - ftp l Vp r a + V &ve 




-2 



-8 V avffl dV r b- r a + 



ir J 




p p 

Dividing through the entire equation by ^roV r puts 

' O 

the equation in the non-dimensional form 
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Momentum Rate ^pj 2 r a 2 ^ ^av^ 2 (r^2-r a 2) _g^ave x ^ / r b“ r a 



Tr«V 2 v2 r 2 v 2 
7T / >Vs o r a Q v s G r a_ v s Q 



V 



,'tr v , 



& I & V.' - 



l£ f r b,- r a j 

2 r r aj 



V 2 

l+T 4 * 



2 2 

.. r 3 






The Momentum Equation is then applied to the ejector in 
the axial direction. Neglecting shear stress at the ejector 
walls, 

Ap rfr„ 2 = ^(Momentum Rate) 



Defining a non-dimensional Momentum Rate by dividing by 

fTfV r % we obtain, 

' © a 0 

Ap r H 2 _ A M 
^ r a 2 ’ 

Defining a non-dimensional pressure difference, Ap*”-, as 

A p» - A P , .All 

¥ V s 2 ‘ 

where Ap=p x -p x=0 and AH *M X 

The final equation used is the Energy Equation in which 
the heat transferred to the surroundings is assumed to be 
zero. The difference in geo-potential energy is negligible 
and, lastly, the change in kinetic energy is negligible in 
comparison with the Enthalpy Rate. Thus the Energy Equation 
becomes 

H -J dll » J" hdm = 2 TTjjCp y* TV rdr =. Constant 
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Using the same profile shapes as before, the Enthalpy 



Rate for a typical axial cross-section becomes 

r r_. p /r-r 
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The second term of the Enthalny Rate Equation, namely, 
r. 



/ bt f v 

r at L 
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Of 



x +dT cos/ r ~ r at 
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-If) rdr 
t L 



was solved by numerical integration. The interval between 

r . and r, ^ was divided into 100 increments and summed, 
at bt 

The other integrals which make up the total value of the 
Enthalpy Rate are integrated and evaluated giving 







In order to place this in non-dimensional form, the 

2 

for numerical integration becomes 



equation was divided by five V r . Thus the integral 
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Non-dimens ionali zing the evaluated portion of the 
Enthalpy Rate gives 
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3. Method 

The equations employed in EJECTMIX I and SJECTMIX II 
were used in the following manner. At a given axial position, 
a core boundary (r a ) is assumed. Primary am’ secondary velo- 
cities (V and V s ) are assumed which will be compatible with 
the Bernoulli Equation for the entrance to the mixing cham- 
ber. A jet boundary (r^) is assumed and the conditions are 
checked for Conservation of Mass by means of the Continuity 
Equation. The primary and secondary velocities are corrected 
by an iterative process until a compatible solution is found. 

The velocities found by the Continuity Equation are then 
used to solve for the pressure at the cross-section under 
consideration. Pressure changes determined from Bernoulli's 
Equation along streamlines in the uniform flow velocity re- 
gions are compared with the pressure change found from the 
overall Momentum Rate. If the pressure is found to be incom- 
patible, the jet boundary (r^) corrected and once again 
the Equation of Continuity is used in order to determine 
compatible central and peripheral velocities. The values 
found by the Continuity Equation are then used to find a new 
value of the Momentum Rate, which, in turn, determines a new 
pressure change which can be used in Bernoulli's Equation to 
check the velocities found by Continuity. This process re- 
peats itself until the core and jet boundaries (r ft and r^), 

the central and peripheral velocities (V and V ) and the 

p s 

pressure change are all compatible in the three bxsic equa- 
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Fig. 4 Block Diagram 



for EJECTMIX I 
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tions used. Then the process is repeated for the next cross- 
section which is displaced axially in the direction of fluid 
flow. 

After several cross-sections have been investigated, the 

core boundary (r a ) or the jet boundary (r^) may reach their 

limiting value, that is, the central core boundary (r ) may 

a 

go to zero or the peripheral jet boundary (r^) may reach the 
mixing chamber boundary (r g ). 

Considering the case of either one of the boundaries 
reaching its limiting value while the other boundary remains 
within its limit, the scheme for finding the compatible solu- 
tion remains exactly the same as that used, when both bound- 
aries were within limits, x-jith the exception of the role of 
Bernoulli's Equation. Since the flow field will no longer 
be uniform x-/here the boundary has reached its limit, the 
streamline along which Bernoulli's Equation Xiras assumed to 
have a constant value will no longer exist. With this in 
mind, Bernoulli's Equation will only be taken into accoxxnt in 
the uniform flow portion which has not reached its limiting 
value. 

When both boundaries have reached their respective limit- 
ing values, no additional restrictions are required other 
than the overall Momentum Rate. Compatible solutions are 
then found using fixed increments of the central velocity as 
the criteria for the particular cross-section, 

EJSCTnIX II is basically the same program as BJ3CTMIX I 
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Fig. 5 Block Diagram for EJECTMIX II 
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with additional restrictions placed upon the flow field. 

These additional restrictions take the form of variations 
in the temperatures of the primary and secondary fluids. The 
temperature boundaries are assumed to spread more rapidly 
than the velocity region boundaries. This has been verified 
by experiments on jets of gases flowing in ejectors.^- Also, 
the temperature is assumed to be a constant until the temp- 
erature boundaries reach their limiting values. The temper- 
ature boundaries are assumed and the value of the energy equa 
tion at the cross-section in question is compared with that 
of the initial cross-section. The central temperature core 
boundary (r is corrected until a compatible solution is 
obtained. 

In order to determine the axial position on which each 

of the representative compatible solutions lie, a spread rate 

E d_ (r^)l is assumed for the peripheral jet boundary. When 
d x j 

this boundary reaches its limiting value, it is assumed that 
the pressure will continue varying as a smooth curve until 
it reaches its limiting value in the mixing chamber. 



^D. Kuchemann and J. Weber, Aerodynamics of Propulsion, 
McGraw-Hill Book Co., Inc., New York, 195>3 
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4. Results 

The data which was obtained using EJECTMIX I and 
EJECTMIX II is presented in graphical form in Appendix 2 and 
in tabular form in Appendix 3» The properties obtained are 
presented in groups according to the ratio of the radius of 
the peripheral boundary to the central jet radius. The 
values of this given ratio are 1.5, 3.0 and 10.0. 

These graphs may be entered with the known physical 
dimensions of the ejector mixing chamber and the velocity 
ratio. The ratio of the initial peripheral radius to the 
central jet radius of ten to one gives an actual area ratio 
of 100 to one and thus approaches the condition whereby the 
peripheral fluid actually would be unbounded. 

The graphical solutions shown are those in which the 
fluid was assumed to return to a uniform flow after it had 
passed doi-m an axial distance equal to ten times the peri- 
pheral diameter. The data and graphs are presented as func- 
tions of axial distance divided by initial central jet radius 

(X/r ). Tabulated data shown in Appendix 3 is the computer 
a 

output which was used to produce these graphs. 

In the course of finding compatible solutions, a temper- 
ature difference of 100 degrees was used. 3y varying only 
the temperature, the solutions obtained for temperature were 
directly proportional to each other. With this in mind, the 
temperature difference for any case can be determined by 
using a proportionality constant and a 100 degree temperature 
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difference. Temperatures are measured above a reference 
temperature which is the secondary fluid temperature at the 
initial cross-section. 

The pressure term as used in Bernoulli's Equation is not 
significant in its absolute value since it cancels out of the 
equation when used in these programs. Its significance lies 
in the change of the non-dimensional pressure term. With 
each specific case an initial non-dimensional pressure can 
be determined. This non-dimensional pressure will be a func- 
tion of the absolute velocities, the pressure and the fluid 
flowing in the ejector. 

Since the initial peripheral velocity ratio is equal to 
one, a change in the non-dimensional pressure of one, before 
the velocity jet boundary reaches its limiting value, leads 
to a zero velocity ratio through the Bernoulli Equation. Any 
increase in the non-dimensional pressure beyond the value of 
one leads to compatible solutions in which the peripheral 
velocity becomes negative. This phenomenon is dependent 
upon the velocity ratio and the cross-sectional area of the 
central and peripheral jets. 

Negative peripheral velocity ratios first appeared in 
the range of the values studied at an initial central velo- 
city ratio of ten and an area ratio less than 100 to one. 

The flow pattern which exists when negative peripheral velo- 
city ratios are generated in the compatible solution is 
shown in Fig. 6. 
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Pig. 6 Profiles and Plow Pattern 
with Negative Peripheral Velocities 



Throughout nost of the program accuracy of one-half of 
one percent is required in order for the solution to be con- 
sidered acceptable. This criterion fails during the transi- 
tion of the peripheral velocity ratio term through its zero 
values. For this portion of the program accuracy of plus or 
minus one one-hundredth for the non-dimensional velocity term 
is required. This change in the criterion ensures that a 
finite and acceptable error can be established for any value 
of the peripheral velocity ratio. 

The final piece of data presented by the program is the 
Momentum Factor. This term gives that multiple of the final 
momentum that exists at any point. 
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5 . Conclusions 

The programs presented here represent the first step in 
a series of programs which would eventually describe the 
physical properties of any fluid flowing in an ejector. This 
logical, continuation would require that work be done on the 
compressible fluid problem at both subsonic and supersonic 
velocities . 

The programs contained in this work solve the ejector 
problem for tvio incompressible fluids of the same density 
whose physical properties are initially different. The data 
obtained is presented here in graphical and tabular form. 
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APPENDIX I 



Appendix I contains a listing, in FORTRAN language, of 
the programs EJECTMIX I and EJECTMIX II. 
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135 VP ( I ) = VP ( I ) + C. CCC1 
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f)P(I) = ( P F T A - .S1GPA - TAL)/(0.5«RS<»2) 
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134 STOP 
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136 R R ( I + 1) = 2.0*RB(1) - RB(I 
GO TO 132 



115 R0( I ) = R 8 ( 1 ) + 0.001 
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128 TAU = VP ( 1 ) ««2«RA( 1 J * 
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APPENDIX II 



Appendix II is a graphical presentation of the data 
obtained by programs EJECTMIX I and EJECTMIX II. Various 
physical characteristics are plotted against the axial dis- 
tance. The axial distance is expressed as a multiple of the 
central jet radius. The various curves shown represent data 
obtained using different initial central velocity ratios. 
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Non-dimensional Secondary Velocity Jet Boundary vs. 
Axial Distance for an Area Ratio of 2. 2£:1 
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Axial Velocity Ratios vs. Axial Distance 
for an Ejector Area Ratio of 2.25:1 
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Peripheral Velocity Ratio vs. Axial Distance 
for an Ejector Area Ratio of 2.25 si 
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Non-dimensional Temperature Jet Boundary 



Non-dimens ional 
Distance for 



Temperature Jet Boundary vs. Axial 
an Ejector Area Ratio of 2.25 ; 1 
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Axial Temperature vs. Axial Distance 
for an Ejector Area Ratio of 2. 2$:1 
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Non-dluionaional Pressure Difference vs. Axial Distance 
for an Ejector Area Ratio of 2.2f>:l 
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Multiple of Pinal Momentum Rate vs. Axial Distance 
for an Ejector Area Ratio of 2.25:1 
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xial Velocity Ratios ( v r _ 0 / v s ) 



Axial Velocity Ratios vs* Axial Distance 
for an Ejector Area Ratio of 9:1 
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Peripheral Velocity Ratio* (V 



Peripheral Velocity Ratio vs. Axial Distance 
; for an Ejector Area Ratio of 9:1 
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Non-dimensional Temperature 
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Axial Temperature vs, 
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Ratio of 9:1 
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Non-dimensional Pressure Difference vs. Axial Distance 
for an Ejector Area Ratio of 9:1 




Axial Distance (X/r a ) 
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Multiple of Pinal Mcnentur. Rate vs 
for an Ejector Area Ratio 



Axial Distance 
of 9:1 




Axial Distance (X/r~) 



66 



Non-dimensional Velocity Jet Boundary 



Non-dimensional Velocity Jet Boundary vs. Axial 
Distance for an Ejector Area Ratio of 100:1 
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Axial Velocity Ratio (V /V 

r -o' 



Axial Velocity Ratio vs. Axial Distance 
for an Ejector Area Ratio of 100:1 




Axial Distance (X/r„) 
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Peripheral Velocity Ratio' (V 



Peripheral Velocity Ratio vs. Axial Distance 
for an Ejector Area Ratio of 100:1 
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Non-dimensional Temperature Jet Boundary 



Non-dimensional Temperature Jet Boundary vs. Axial 
Distance. for an Ejector Area Ratio of 100:1 
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Axial Temperature (T 



Axial Temperature vs* Axial Distance 
for an Ejector Area Ratio of 100:1, 
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Non-dimensional Pressure Difference vs. Axial Distance 
for an Ejector Area Ratio of 100:1 
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Multiple of Final Momentum Rato, 



Multiple of Final Momentum Rate vs. Axial Distance 
for an Ejector Area Ratio of 100:1 




APPENDIX III 



Appendix III contains edited data from which the graphs 
of Appendix II were drawn. The headings are the sane as 
defined in the section of Symbols and Abbreviations, 
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Edited Data Tor Initial Conditions of 
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Edited Data for Initial Conditions of 
Central Velocity Ratio of 1,5 and Peripheral Radius Ratio of 3»0 
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Edited Data for Initial Conditions of 
Central Velocity Ratio of 3.0 and Peripheral Radius Ratio of 3.0 
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